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Abstract 



We propose a deterministic weighted scale-free small-world model for considering 

■ pseudofractal web with the coevolution of topology and weight. Considering the 
^ . fluctuations in traffic flow constitute a main reason for congestion of packet delivery 
Q I and poor performance of communication networks, we suggest a recursive algorithm 

to generate the network, which restricts the traffic fluctuations on it effectively 

■ during the evolutionary process, we provide a relatively complete view of topological 
[ structure and weight dynamics characteristics of the networks: weight and strength 

j. I distribution; degree correlations; average clustering coefficient and degree-cluster 

correlations; as well as the diameter. 

Key words: Complex networks. Scale- free networks. Weighted networks, 
Q\ [ Disordered systems. Traffic fluctuations 
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1 Introduction 



To understand the general principles in architectures of networks, many deter- 
ministic models are introduced into complex networks 

* Corresponding author. 

Email addresses: zhangzz@fudan.edu.cn (Zhongzhi Zhang), 
sgzhou@fudan.edu.cn (Shuigeng Zhou), jhguan@mail.tongji.edu.cn (Jihong 
Guan). 



Preprint submitted to Elsevier 



These models are useful tools for investigating analytically not only topological 
features of networks in detail p!|2|3|^f5|6|7f8|9|10fllfl2] . but also dynamical 
problems on the networks [T^14fl5|ITB] . Before presenting our own findings, it 
is worth reviewing some of this preceding work to understand its achievements 
and shortcomings. Deterministic scale-free networks were firstly proposed by 
Barabasi et al. in Ref. [1] and intensively studied in Ref. |2] to generate a 
scale- free topology. However, to some extent, the small exponent 7 of the de- 
gree distribution for the model didn't satisfy the real statistic results well. 
Instead, Dorogovtsev et al. introduced another elegant model, called pseud- 
ofractal scale-free web (PSW) |3] which is extended by Comellas et al. conse- 
quently [1]. Based on a similar idea of PSW, Jung et al. presented a class of 
recursive trees [5] , which have the small- world behavior built in. Additionally, 
in order to discuss modularity, Ravasz et al. proposed a hierarchical network 
model pi7] . the exact scaling properties and extensive study of which were re- 
ported in Refs. [S] and [2], respectively. Recently, motivated by the problem of 
Apollonian space-filing packing, Apollonian networks [10] with a typical loop 
structure were introduced and intensively investigated [T7|18|19|20|21] . These 
pioneering works are all invaluable tools for the topology of networks studies. 

In the last few years, it is found that many real networks are inhomogeneous, 
consisting of distinct nodes and links. For instance, the scientist collabora- 
tion network, where scientists are identified with nodes, and an edge exists 
between two scientists if they have coauthored at least one paper [22], and the 
Internet at the AS level, where the link weights represent the bandwidth of a 
cable and node weight the load of a router [23], among other areas. Recently, 
weight dynamics ideas have been applied with success to topics as diverse, such 
as random walks [21], condensation [25], synchronization [26], traffic conges- 
tion [27], epidemic spreading [28|29] . information filtering [30], to name but a 
few. The findings above might provide insight for understanding the correla- 
tions among weighted quantities and the underlying topological structure and 
dynamics behaviors of the weighted networks. 

Most previous weighted random models [3T|32f33ll34f35] with topology and 
weight coevolution, however, possess very loose clustering structures when the 
size of the networks is large. At the same time, previous deterministic models, 

which ignore the heterogeneity of edges in real networks. What's more, the 
models [3lH] on PSW networks fail to provide the reason for adopting the recur- 
sive way to build up the networks. Consequently, in this paper, we introduce 
a model bringing weight evolution into the growth of pseudofractal scale-free 
web (PSW) [3] that aims to circumvent these incongruities properly. As we 
will show, in the case of the recursive construction, the traffic and its fiuctua- 
tion decrease exponentially with time either on edges or at nodes. Hence, we 
believe the construction method may shed some light on networks design to 
improve the control and speed of the whole network [SB]. At the same time. 
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our comprehensive and rigorous solutions may help people understand better 
the interplay between network topology and weight dynamics. 



2 The model 

The construction of the model is controlled by two parameters m and 5, evolv- 
ing in a recursive way. We denote the network after t steps by G(t), t > 
(see Fig. [1]). Then the network at step t is constructed as follows. For t = 0, 
R{0) is a triangle consisting of three links with unit weight. For t > 1, G{t) 
is obtained from G(t — 1). We add mw (m is positive integer) new nodes for 
each of the links with weight w, and we connect each new node to both ends 
of this link by new links of unit weight; moreover, we increase the weight of 
these links by m5w {6 is positive integer). 

Before introducing our model further, we explain why adopting such a recur- 
sive way and why the generated network is increasingly efficient for transmit- 
ting information with network order. In this model, the recursive construction 
is motivated by the practical need to improve the transport capacity of real 
networks. As is known to us, both the physical networks and the numbers of 
users are growing continuously. The performance of the networks for larger 
system sizes and heavier loads are critical issues to be addressed in order to 
guarantee networks' functioning in he near future. For example, if the traffic 
fluctuates dramatically, a highway is more likely to be congested frequently 
when the peak value of traffic exceeding its capacity. There is thus a need to 
build up more branches to distribute the heavy traffic. However, "how many 
branches must we have?" and 'Where shall we put them?" are open questions 
yet. 

Recently, the authors of Ref. [37] claimed the fluctuations in traffic flow consti- 
tute the main factor affecting the performance of networks. They derived the 
dependence of fluctuations with the mean traffic on unweighted networks ana- 
lytically. Consequently, their recipes were adopted extensively to the weighted 
networks by the authors of Ref [3B] • As shown in Ref [5H] , for correlated net- 
works (assortative or disassortative mixing [3310]), the average traffic through 
a link Lij during a time window can be represented as 




(1) 



and the standard deviation can be expressed as 
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Fig. 1. Illustration of the deterministically growing network for the particular case 
of m = 2 and 6 = 1, showing the first three steps of growing process. The gray links 
in the figure denote the links with weight 1, the red links with weight 3, and the 
blue links with weight 9. The number of rings around a gray node denotes its age. 

where Wij is its link weight [3H]. The length of the time window for observation 
is M. The average number of cars or walkers among various time windows in 
the network is denoted by R. A is defined as a random variable representing 
the number of walkers travelling through the link in the time window. 

With respect to the traffic at nodes, the average traffic at a node i 

if.) = ^^RM, (3) 

Then the standard deviation as a function of (/j) can be obtained as: 

+ (4) 

At each time step, the traffic will be distributed dispersedly to the newly built 
edges (nodes). The larger, the size of the considered network is, the lower the 
fluctuations on each edge (at each node) should be. Details of the analysis are 
provided in section 13.11 and 13.21 

Notice that there are in fact three limiting cases of the present model. In the 
special case m = 1 and 5 = 0, it is reduced to the pseudofractal scale-free web 
described in [3]. When 5 = 0, it is a particular case of the geometric growth 
networks discussed in [H]. When m = 1, it is the same as the deterministic 
weighted networks proposed in [12]. Thus, vary parameters m and 6, we can 
study many crossovers between these limiting cases. 

Let us consider the total number of nodes N^, the total number of links 
and the total weight of all links Wt in G{t). The number of nodes created at 
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step t is denoted by n„(t). Note that the addition of each new node leads to 
two new hnks, so the number of hnks generated at step t is ne(t) = 2n^{t). 
By construction, for t > 1, we have 

= mWt_i, (5) 



Et = Et_i + 2n,{t), 



(6) 



and 



Wt = Wt-i{l + m6) + 2mWt-i. (7) 

On the right-hand side of Eq. ([7]), the first item is the sum of weight of the 
old links, and the second term describes the total weight of the new links 
generated in step t. Eq. ([7]) can be simplified to 



Wt = {l + m6 + 2m)Wt-i. 
Considering the initial condition Wq = 3, we obtain 

Wt = 3{l + m6 + 2m)*. 



(8) 



(9) 



Substituting Eq. (|9]) into Eq. (|5]) and using Wq = 3, the number of nodes 
created at step t (t > 1) is obtained to be 



ny{t) = 3m(l + mS + 2m) 



t-1 



(10) 



Hence, one can figure out the growth of the network is accelerated. Then the 
total number of nodes present at step t is 



ti=0 



2 + S 



[l + mS + 2m)* + S + 1 



Combining Eq. ( fTOj) with Eq. ([6]) and considering i?o = 3, it follows that 



Et 



2 + 6 



'l + mS + 2m)* 



6 



;i2) 



Thus for large t, the average degree kt = ^ is approximately equal to 4. 



3 Structural properties 



In what follows we will study how the tunable parameters m and S control 
some relevant characteristics of the weighted network G{t). Firstly, we give 
out the analytical solution of distribution of strength, degree and weight to 
test its scale-free nature; simultaneously, we show the analytical expression of 
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average traffic and its deviation of nodes and edges; subsequently, we move 
forward to the average clustering coefficient coupled with the diameter of this 
network for the sake of verifing its small-world property; finally, we study the 
degree correlations as well. 



3.1 Weight distribution 



Let We{t) be the weight of link e at step t. In the view of all the links emerging 
simultaneously have the same weight, it can be recast recursivel as follows 



w^{t) = {l + m6)we{t - 1). 



(13) 



If link e enters the network at step r, then Weir) = 1. Thus, we can easily 
have 

We{t) = {1 + m6y-\ (14) 

Obviously, the weight spectrum of the network is discrete. It follows that the 
weight distribution is given by 



P(w) 



2 + 6 



2{1 + mS + 2my + 



Heir) 2m{2 + 6) {1 + m6 + 2m) 



2(1 + mS + 2mf + 







for r = 0, 

for r > 1, (15) 

otherwise 



and that the cumulative weight distribution [43)13] is 

r2e(/i) 2(l + m(5 + 2m)^ 



M = E 



E, 



2(1 + m(5 + 2m)* + | 



Substituting for r in this expression using r = t — gives 



(16) 



ln(l + mc5 + 2m) 

2(1 + m(5 + 2m)*w in(i+ma) +| 
^^"'"^'^^ ^ 2(l + m5 + 2m)* + | 

ln(l+mi + 2m) 

in{i+ms) for large t. (17) 



Apparently, the weight distribution follows a power law with the exponent 
7^ = 1 + ^^^ij^l^^s)"^ ■ For the particular case of m = 1, Eq. ( [T7|) recovers the 
result previously obtained in Ref. |12] 



In this paper, the networks generated by the model are disassortative, which 
will be discussed analytically in section 13.61 For disassortative networks, the 
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fluctuation aij on an edge Lij relies on average traffic fij, while fij is governed 
by with R and M fixed. Inserting Eq. ([14]) and Eq. into Eq. (II]), 

one can obtain the average traffic on an arbitrary edge Lij can be written as 

^^'^^ - ^WT " 3(l + m5 + 2m)*^^' ^^^^ 

where Tjj denotes the entry time of Lij. Thus the standard deviation can be 
expressed as 
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l+m5 V {I + m5)-^'i RM 



1 + m5 + 2m y 3 

l + m5 V* (1 + m(5)-2-»^ (A^ + A)M2 



1 + m(5 + 2m / 27 



(19) 



Eq. f|T8l) and Eq. f|T9]) show that both and (Tjj decrease as an exponential 
function of time t with R and W fixed as a result of ^^T",^^ < 1. Notice that, 

l+md+2m ' 

we consider the parameter i? as a const with the time evolution in our model 
in that the relation between R{t) and t depends on various specific systems, 
which is not the focus of the present paper. In fact, one can easily observe that 
the recursive algorithm can also restrict the traffic fiuctuations on edges in the 
considered networks effectively during the evolutionary process, in which R{t) 
is not larger than C~^^^g"^ y- Discussion of traffic and its fluctuation at nodes 
will subsequently be given in the next section. 



3.2 Strength distribution 



In a weighted network, a node strength is a natural generalization of its con- 
nectivity. The strength Si of node i is defined as 

Si= ^ij ' (20) 

where Wy denotes the weight of the link between nodes i and j, Qi is the set 
of all the nearest neighbors of i. The strength distribution P{s) measures the 
probability that a randomly selected node has exactly strength s. 

Let Si(t) be the strength of node i at step t. If node i is added to the network 
at step ti, then Sj(tj) = 2. Moreover, we introduce the quantity Asi(t), which 
is defined as the difference between Si(t) and Sj(t — 1). By construction, one 
can easily obtain 
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Asi(t) = Si{t) - Si{t - 1) = mS Wij + m ^ 



= m6si{t-l) + msi{t-l). (21) 

Here the first item accounts for tlie increase of weight of the old hnks existing 
in step t — 1. The second term describes the total weigh of the new hnks with 
unit weight that are generated in step t and connected to i. 

From Eq. (12T|) . one can derive following recursive relation: 

Si{t) = {l + m5 + m)si{t-l). (22) 

Using Si(ti) = 2, we obtain 

Si(t) = 2(1 +m5 + m)*-*\ (23) 

Since the strength of each node has been obtained explicitly in Eq. ( I23l) . we 
can get the strength distribution via its cumulative distribution [3P3] . i.e. 

p / N _ V- M/f) _ jl + mS + 2mY^ + S + 1 
ku ^* - {l + m6 + 2my + 5 + l- 

From Eq. (|23|) . we can derive ti = t — ij^(l^^s+m) ■ Substituting the obtained 
result of ti into Eq. 1^^, we have 



ln(l + mg + 2m) 

(l + m(5 + 2m)* (f) '"(^+™^+'") +6 + 1 
{1 + m6 + 2my + 6 + 1 

ln(l + m(S + 2m) 
S \ ~ ln(l + m6+m) 

-1 for large t. (25) 

Thus, node strength distribution exhibits a power law behavior with the ex- 
ponent 7^ = 1 + ^^^^i^s+m) • ^'-'^ special case m = 1, Eq. (125!) recovers the 
results previously reported in Ref. [12]. 



On the other hand, the fluctuation at an arbitrary node i is based on ^^^^^ — . 

Inserting Eq. fl23|) and Eq. (|9]) into Eq. ([1]), one can obtain the average traffic 
at an arbitrary node i can be represented as 

2Wt 3{l + m6 + 2my ' ^ ^ 

where ti denotes the entry time of the node i. Then the standard deviation as 
a function of (/,) is 
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mS + m \ * (1 + mS + m) *'RM 



mS + 2m I 3 

2t /- 



m5 + m)-^*'(A2 + A)M2 ^^^^ 



mS + 2m I 27 



It is easy to find that both /j and cTj decrease exponentially with t as -|^'^^/_^2m ^ 
1, which is similar with the former result on edges. Although the strength 
of nodes growing exponentially, the traffic and fiuctuation at them still de- 
crease with the growing size of the networks in the case that i? is a con- 
stant or R{t) < C'^^^g^^ Y- In other words, the sufficient condition of keeping 
the potential traffic fluctuation problems away from the resulting networks is 

R{t) < C'^^j^s"^ ^ ir,ii+ms+2m) qy the average number of walkers is invariable. 
The novel property is interesting and has not been investigated by previous 
works [3|31f32f33f34f35] . Therefore, to some extent, this model may provide a 
paradigm to control the traffic fluctuations and improve transport efficiency 
of the whole network 



3.3 Degree distribution 



Similarly to the strength, all simultaneously emerging nodes have the same 
degree. Let ki{t) be the degree of node i at step t. If node i is added to the 
network at step ti, then by construction fcj(tj) = 2. After that, the degree ki{t) 
evolves as 

ki{t) = ki{t-l) +ms,{t-l), (28) 
where msi{t—l) is the degree increment Aki{t) of node i at step t. Substituting 
Eq. ([23]) into Eq. we have 



Aki{t) = 2m (1 + mS + my-^-^\ (29) 
Consequently, the degree ki{t) of node i at time t is 

* 2 r 1 

ki{t) = ki{t,) + ^Hv) = '2 + TVT (mS + 1 + my-'^ - 1 . (30) 

Analogously to computation of cumulative strength distribution, one can find 
the cumulative degree distribution 



-Pcum(^) 



{l + m6 + 2my f (5 + 1) - 5 



ln(l + m(5 + 2m) 
ln{l + m(5 + m) ^ 



{l + m5 + 2m)* + 6+1 



k 



(6 + 1) -6 



ln(l + m<5 + 2772) 
ln(l + m(5 + m) 



for large t. 



(31) 
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Fig. 2. Cumulative degree distribution Pcum{k) versus k for different 6 and m cor- 
responding to Eq. (j3ip . The measurements are taken at t = 5, illustrating these 
networks display a power-law degree distribution. The dashed lines are the best fits, 
with -fcum = 1.20806, -1.12887, -1.15699 respectively. 

As is shown in the Fig. [21 the degree distribution has the scale-free property 
with the same exponent as 7., {'jk = 7s = 1 + 7c«m, where 'jcum = ^^TT^^iB^) ■ 



3.4 Clustering coefficient 



In this model, the analytical expression for clustering coefficient C{k) of the 
individual node with degree k can be derived exactly. For instance, if a new 
node is connected to both ends of a link, its degree and clustering coefficient 
will be 2 and 1, respectively. Naturally, its degree will increase by one when 
connecting a new node in the next step. On the other hand, there must be 
an existing neighbor of it attaching to the new node at the same time. On 
the other hand, there must be an existing neighbor of it, attaching to the new 
node as well. Because our networks are corresponding to the particular case 
g = 2 of the recursive clique trees [3114] . for a node of degree k, we have 

C{k) = = 2/k. (32) 

2 

The scaling C( A;) ~ k^^ has been found for some network models [71ll0fl2fl7lHl] . 
and has also observed in several real-life networks [7]. 
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Fig. 3. The solutions of Eq. ()33p for 6 and m ranging from 1 to 5 respectively. The 
measurements are taken att = 100, illustrating these networks display a high degree 
of clustering. 



Using Eq. (132|) . we can obtain the clustering Ct of the networks at step t: 

^-.-^i:^. (33) 

where the sum is the total of clustering coefficient for all nodes and kr = 
2 + [{m6 + 1 + my~^ — 1] shown by Eq. ( 15U]) is the degree of the nodes 
created at step r. 

It can be easily proved that Ct increases with q for arbitrary fixed m, and likely 
Ct increases with m when q fixed. In the case of t = 100 (A^ — )■ oo), Eq. (l33l) 
converges to a correspondingly large value C. When 5 = 2, for m = 1,2,3,4 
and 5, C equal to 0.886,_0.922, 0.941, 0.952 and 0.96. respectively. When m = 2, 
for (5 = 1, 2, 3, 4 and 5, C are 0.899, 0.922, 0.937, 0.947 and 0.954, respectively. 
Evidently, the clustering coefficient of our networks is correspondingly stable 
and close to 1. Moreover, the average clustering coefficient C can be tuned by 
S and m (see Fig. [3]) . 

In the classical weighted co-evolutionary models, for example, BBV networks 
[22|32] . the average clustering coefficient rapidly decreases when the networks 
is growing (see Fig. H]). Our simulations confirm that in the limit of large 
networks [N ^ 1), the BBV networks' clustering coefficient is getting close to 
zero. However, many real-world networks have a relatively stable and nonzero 
clustering coefficient, which make the results of BBV model useful, but far 
from comprehensive. In the Fig. |5l we performed numerical solutions for our 
model at various values of S with fixed m in the panel (a) (various values of m 
with fixed 6 in the panel (b)) up to t = 50. For the infinite network, one can 
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Fig. 4. Average clustering coefficient C as a function of size N of nodes for different 
6 in the BBV networks are shown in the inset, using the weight of a new hnk wq = 1, 
the size of initial seed niQ = 3, the degree of a new node m = 3. 
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Fig. 5. Semilogarithmic plots of the average clustering coefficient C against net- 
works size N for (a) 6 = 1,2, ...,6 when m = 2 and (b) m = 1,2, ...,6 when 
6 = 2 corresponding to Eq. (j33p . The measurements are taken over a time interval 
t = 1,2, 50, illustrating that C is a finite nonzero value independent of N for our 
model. 

easily obtain that this tunable average clustering coefficients of our networks 
is independent of their size and tends to a nonzero limit, which is a unique 
property shared by many real networks pH3] . At the same time, this feature 
gives excellent agreement with the previous unweighted studies [3]. 



3.5 Diameter 



As is known to all, the diameter of a network is defined as the maximum of the 
shortest distances between all pairs of nodes, which characterizes the longest 
communication delay in the network. Small diameter is an important property 
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of small-world network. Fortunately, for our networks, it can be calculated 
easily. Below we give the precise analytical computation of diameter of G{t) 
denoted by Diam{G{t)). 

It is easy to see that at step t = (resp. t = 1), the diameter is equal to 
1 (resp. 2). At each step t > 2, one can easily see that the diameter always 
lies on a pair of nodes created at step t. In order to simplify the analysis, 
we first note that it is unnecessary to check all the nodes in the networks 
to fix the diameter. In other words, some nodes ("inner" nodes) added at a 
given step can be ignored, because they do not increase the diameter of the 
previous net. Here, so-called "inner" nodes are those that connect to links 
that already existed before step t — 1. Indeed, for these nodes we know that a 
similar construction has been done in previous steps, so they have nothing to 
do with the calculation of the diameter. 

Let us call "outer" nodes the nodes which are connected to a fresh link. Clearly, 
at each step, the diameter depends on the distances between outer nodes. At 
any step t > 2, we note that an outer node cannot be linked with two nodes 
created during the same step r < t — 1. Indeed, we know that from step 2, 
no outer node is connected to two nodes of the initial triangle G{0). Thus, 
for any step t > 2, any outer node is connected with nodes that appeared at 
pairwise different steps. Now consider two outer nodes created at step t > 2, 
say Vt and Wt- Then Vt is connected to two nodes, and one of them must have 
been created before or during step t — 2. 

We summarize the above arguments, and gather them into two cases: (a) t = 21 
is even. Then, if we make / 'jumps' from Vt we reach the initial triangle G'(O) 
in which we can reach any Wt by using a link of G{0) and making / jumps 
to Wt in a similar way. Thus Diam(G{2l)) < 2/ + 1 = t + 1. (b) t = 2/ + 1 
is odd. In this case we can stop after / jumps at G'(l), for which we know 
that the diameter is 2, and make / jumps in a similar way to reach Wt- Thus 
Dia'm{G{2l + 1)) < 2{l + 1) = t + 1. Obviously, the bound can be reached by 
pairs of outer nodes created at step t. More precisely, these two nodes Vt and 
Wt share the property, that both of them are connected to two nodes added 
at steps t — 1, t — 2 respectively. Hence, formally, Diam{G{t)) = t + 1 for any 
t > 0. Considering Nt ^ {1 + m6 + 2m)*, the diameter is small and scales 
logarithmically with the number of network nodes. 



3.6 Degree correlations 

In complex network, degree correlations [39]|40f45f46ll47j . has attracted much 
attention, because it can give out a unique description of network struc- 
tures, which could help researchers understand the characteristics of net- 
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works [39|46f47|H8] ■ An interesting quantity related to degree correlations is 
the average degree of the nearest neighbors for nodes with degree k, denoted 
as knn{k) [1^471139] ■ When fc„„(/c) increases with fc, it means that nodes have 
a tendency to connect to nodes with a similar or larger degree. In this case 
the network is defined as assortative [39lH0] . In contrast, if knn{k) is decreas- 
ing with k, which implies that nodes of large degree are likely to have near 
neighbors with small degree, then the network is said to be disassortative. If 
correlations are absent, knn{k) = const. 

In our networks, we can acquire k^nik) exactly by Eq. (1501) . Except for three 
initial nodes generated at step 0, no nodes born in the same step, will be linked 
to each other. All links from the newcomers to old nodes with larger degree 
are made at their creation steps. Then, these newcomers become old ones to 
accept the nodes with smaller degree made at each subsequent steps. These 
results are shown in the expression 



knn{,k) = ^ A rn-n^{t'i)s{t[ji-l)k{t'i,t) 

tr=t 

+ m-n,{ti)s{U,t',-l)k{t',,t)). (34) 

i-=ti+l 

Here the first sum on the right-hand side accounts for the links made to nodes 
with larger degree (i.e., t[ < ti) when the node was generated at ti. The second 
sum describes the links made to the current smallest degree nodes at each step 
t', > U. 

Substituting Eqs. (fTOj) and ( 130|1 into Eq. ( 134|) . one expects that 



, {m5 + I + 2m){m5 + I + m) 

mo^ + 2mo + m + o 
(m5 + l+m)2*' t-ti 

+ r , 1 , ]' (35) 



(m6 + 1 + 2m)** m5 + 1 + m 

in the infinite limit of t, where kr ~ s^i^^ ~^ 1 -|-m)*~*\ In another word, the 
initial step knniti,t) grows linearly with time. Consequently, writing Eq. ( j35l) 
in terms of k, it is straightforward to obtain 



knn{k, t) 



^ {m5 + 1 + 2m) {m5 + 1 + m) {m6 + 1 + m)^* 
mS'^ + 2m6 + m + 6 {m5 + 1 + 2m)* 



k-2 



{5 + l) + l 



2 ln(l + m(5 + m)-ln(l + m6 + m) 
ln(l + m(5 + m) 



(36) 
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Apparently, knn{k) is approximately a power law function of k with nega- 
tive exponent, which indicates that the networks are disassortative. Note that 
knn{k) of the Internet exhibits a similar power-law dependence on the degree 
knnik) ~ k-"", with w = 0.5 HE]. 



4 Conclusion and discussion 

To sum up, we have proposed and investigated a deterministic weighted net- 
work model, which is constructed in a recursive fashion. The recursive con- 
struction guarantees that the traffic fluctuations of nodes and edges decrease 
exponentially with the time of evolution. The weights of these networks char- 
acterizing the various connections exhibit complex statistical features with 
highly tunable degree, strength, and weight distributions, which display power- 
law behavior. We have shown the analytical results for degree distributions 
with tunable exponent and large clustering coefficient, as well as small diam- 
eter. Particularly, the features of clustering coefficient in our proposed model, 
i.e., it is independent of its net size, might lead to a better understanding of 
realistic networks. To some extent, our model can thus perform well in con- 
trolling and designing a variety of weighted scale-free small-world networks to 
improve their transport efficiency. 
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